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We study a model of strongly correlated S — 1/2 fermions on the planar pyrochlore, or checker- 
board, lattice, at fractional (1/8) filling. Starting with the extended Hubbard model in the limit 
of strong Coulomb repulsion, low-energy configurations can be mapped onto hard-core dimer con- 
figurations whose dimers carry a spin degree of freedom. An effective Hamiltonian similar to the 
kinetic term of a quantum dimer model on the square lattice which rotates two parallel dimers (in a 
spin-singlet configuration) by 90°-degrees naturally emerges. We also introduce an additional term 
in the Hamiltonian, a generalized dimer plaquette interaction, in order to realize a closer analogy 
to the latter model. For a strong dimer plaquette attraction stabilizing a columnar phase, a spon- 
taneous dimerization takes place in the direction of the columns of (spin-carrying) dimers. Using 
exact diagonalizations of two-dimensional periodic clusters, the analysis of the low-energy spectrum 
and of several types of correlation functions gives indeed evidence for a new type of lattice symmetry 
breaking phase, the eight-fold degenerate Mixed Columnar-Plaquette Crystal, and for a transition 
from this phase to a Resonating-Singlet-Pair Crystal (found in previous studies) which restores the 
rotational symmetry of the lattice. Similar conclusions and phase diagram are also reached from a 
simple variational approach. 



I. MODEL, PURPOSES, AND METHOD 

A. Introduction and summary of previous results 

The interplay between electronic correlations and the 
lattice geometry in quantum magnets can lead to a rich 
variety of spin gapped disordered phases, either spin 
liquids with fractionalized excitations or various types 
of Valence Bond Crystals (VBC), which break sponta- 
neously some of the lattice symmetries. Among ma- 
terials magnetically frustrated and possibly presenting 
such phases, those with a pyrochlore structure, a three- 
dimensional (3D) array of corner-sharing tetrahedra, are 
of particular interest because of the absence of mag- 
netic order down to very low temperatures [lj. On 
a two-dimensional (2D) version of the pyrochlore lat- 
tice, the checkerboard lattice (see Fig.QJ, the Hciscnberg 
model presents a VBC of particular interest, the plaque- 
tte phase Q, which exhibits the rotational symmetry of 
the lattice. To understand the physics of undoped and 
doped frustrated magnets and predict the occurrence of 
these phases in real materials, theoretical tools such as 
the Hubbard model, and models derived from it in the 
limit of strong on-site repulsion, are commonly used in 
2D (also in 3D) systems. These exhibit very interest- 
ing properties: in a model of bosons on the triangular- 
lattice, doping away from commensurate fillings drives a 
transition from an insulator to a supersolid phase (with 
charge ordering and a finite superfluid density) [|| , which 
is also found in a model of bosons on the checkerboard 
lattice d|; on the same lattice, spinless fermions near 
1 /4 filling present interesting properties such as fractional 
charge excitations 0, 0| . 



In a more specific context, to describe the non- 
magnetic Resonating Valence-Bond phase of cuprate ma- 
terials and the transition of this phase to the super- 
conducting phase, the Quantum Dimer Model (QDM) 
was developed in the late eighties Q, mostly on two- 
dimensional lattices. This model displays different types 
of Valence Bond Crystals (with close analogies with their 
spin counterparts), among which the quite exotic pla- 
quette phases, and, depending on the lattice (non-) bi- 
partitcness, cither a liquid phase with topological order 
(on the triangular lattice @ ) or a quantum critical point 
(Roksar-Kivelson point), both presenting deconfined ex- 
citations Q. The QDM is also connected to the physics 
of pyrochlore systems, since strong Coulomb repulsion 
in the (extended) Hubbard model on the kagome or on 
the checkerboard lattice, either for bosons or fermions at 
special fractional fillings, select low-energy configurations 
that can be mapped onto dimer or loop configurations. 

The present study belongs to a series of works about 
the extended Hubbard model on the 2D pyrochlore (or 
checkerboard) lattice, at and slightly away from frac- 
tional fillings [!, H, [l(| [lH ■ Here we focus on the effective 
model of S = 1/2 fermions at 1/8 filling (1 particle for 4 
sites) derived from the extended Hubbard model, 
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in the limit of large Coulomb repulsions U and V (with 
respect to the energy scale given by the particle hopping 
amplitude t). Let us briefly review our current under- 
standing of the physics of the extended Hubbard model 
on the pyrochlore lattice at the special fractional fillings 
n = k/4 (k = 1,2,3) 0, LT3] and mention the important 
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FIG. 1: Schematic representation of the Hubbard model on 
the checkerboard lattice. In the limit \t\ <C V <C U, this model 
becomes equivalent to that described by the Hamiltonian Hk 
(see Eq. \T% 



remaining issues. For spinful fermions at these fillings, 
in the limit where the on-site repulsion U is very large 
compared to the nearest neighbour repulsion V and hop- 
ping £, a mctal-to-insulator transition was found for in- 
creasing \V/t\ at the filling factor n = 1/4 (similarly as 
in the corresponding hard-core bosonic model (HI), an d 
the corresponding insulating phase (for |£| -C V -C U) ex- 
hibits plaquette correlations indicating an ordering very 
different from a simple charge density wave. Whether 
the metal-to-insulator transition occurs immediately at 
infinitesimal V or at a finite value depends on the de- 
gree of the commcnsurability k. Indeed, the perfect nest- 
ing property of the non-interaction Fermi surface realized 
only for k = 2 (and for a given sign of £) leads to an in- 
stability for arbitrarily small U and V In the limit 
of interest here (strong couplings), an effective model of 
S = 1/2 fermions was derived in the same study, in- 
volving a 2-particlc hopping term (amplitude £2) and an 
additional term (amplitude W) counting the number of 
singlet pairs on uncrossed plaquettes. Varying the ra- 
tio of the amplitudes of these terms, the system can be 
tuned from a charge ordered columnar phase (the inter- 
nal structure of columns being that of Hciscnbcrg AF 
chains) in the limit W -C — 1*2 1 to a disordered RK point 
at W = ti- Note that, in this formulation, the case 
W = is believed to provide the effective description 
of the insulating phase of the largc-U, large- V Hubbard 
model on the checkerboard lattice mentioned above. 

So far the case corresponding to filling n = 1/2 is un- 
derstood the best. A phase transition was clearly ev- 
idenced at finite (negative) W '/t% between the charge- 
ordered phase and a Resonating Singlet-Pair Crystal 
(RSPC) using an analysis based on the symmetry- 
resolved low-energy spectrum and plaquette correlations 
in the ground state (GS). The system at n = 3/4 was 
shown to exhibit also plaquette order by forming a (lat- 
tice rotationally-invariant) Resonating Singlet-Pair Crys- 
tal, although with a quadrupling of the lattice unit cell 
(instead of a doubling for n = 1/2) and a 4- fold degener- 
ate ground state. Concerning the n = 1/4 case, the con- 
clusions of the previous study were less clear: although 



the evolution of plaquette correlations with W ' jti also 
supports a transition from a charge-ordered to a RSPC, 
the analysis of low-energy eigenstates was less conclusive 
than in the n = 1/2 and n = 3/4 cases, primarily due to 
larger finite-size effects: in the previous study computa- 
tions were done on a N = 32 checkerboard cluster with 
periodic boundary conditions. Moreover for n = 1/4, 
taking into account the possibility of new mixed phases 
(which are not charge-localized but break rotational sym- 
metry) requires more caution in the analysis of the low- 
energy spectrum (and hence larger clusters). This leads 
us to consider a new scenario for the phase diagram of the 
model, which will be described in more detail hereafter. 

B. The effective model 

As outlined above, the effective model is derived from 
the extended-Hubbard model for S = 1/2 fcrmionic par- 
ticles on the checkerboard lattice, in the limit of very 
large Coulomb repulsion (more precisely |£| <C V <C U). 
In this limit, at 1/8 filling (n = 1/4), one can exclude 
configurations where the 2 neighboring sites arc simulta- 
neously occupied. In other words, each tctrahedra should 
contain exactly one particle (of cither spin), an ice rule- 
type constraint, which still leaves an exponentially large 
number of states. As discussed in the literature, once the 
particles are viewed as dimers living on the bonds of the 
square lattice formed by the centers of crossed plaque- 
ttes, this constraint is equivalent to the hard core dimcr 
constraint on the square lattice. However, in contrast to 
the "usual" QDM on the square lattice, here each dimcr 
carries a color index (associated to the spin of the elec- 
tron it represents). 

In this limit, a single particle hops out of a low-energy 
configuration (colored dimer configuration) creates a de- 
fect centered on a tctrahedra with an energy cost V . This 
defect can however be annealed by the subsequent hop- 
ping of the second particle on the " defect tetrahedron" . 
Such processes lead to an effective kinetic term, i.e. a 
correlated 2-particle hopping, of amplitude £2 = 2t 2 /V. 
In terms of dimers, this term looks like the kinetic term 
of the RK model, but acts only on particles of opposite 
spin on the same uncrossed plaquette (i.e. dimers of op- 
posite color on the same plaquette). The particles being 
fermionic, the expression of the kinetic term involves op- 
erators of creation (destruction) of singlets on uncrossed 
plaquettes, c^cj^ - c^c] T (c tT c,x - Cjjc^), 

H K = -t 2 ( ( C 1 T C ]| - 4i c tf) x 

<ijkl> \ 

(c fcT c u - c kl c n ) + c.c^j , (1.2) 

where the sum is on uncrossed plaquettes (going around 
a plaquette < ijkl >, sites are in the order i,k,j,l). A 
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unitary transformation, consisting in defining operators 
&jt = — c-j on every other ascending and every other 
descending line of the checkerboard lattice oriented as 
on Fig. [T] (i.e. every other vertical line of vertical links 
and every other horizontal line of horizontal links of the 
associated square dimer lattice) and &vj = cfl otherwise, 
allows that each 2-particle hopping term to have the same 
amplitude — t 2 in terms of b^l operators, 

<ijkl> \ 

(bjfct&u + b kl b^) + c.c^j . (1.3) 

Notice that this is valid only in the insulating phases at 
specific fractional fillings like n — 1/4, thanks to the ice 
rule-type constraint. In addition, it is possible to label 
the sites of the lattice in such a way that all exchange 
processes on the empty squares do not involve any re- 
ordering of the fermions so that the b^l operators can be 
considered as bosonic. In other words, our new formula- 
tion uses the bosonic representation of the spin singlets 
(we have checked the equivalence numerically on the 32- 
site cluster). 

Following the initial suggestion of Ref. Q and accord- 
ing to the discussion above, we also consider a term anal- 
ogous to the potential term of the QDM, although here 
it is no longer diagonal in the basis of configurations, 

H w = W J2 [n i n j (l-n k ){l-ni)(l/2-2S i .S j ) 

<ijkl> \ 

+ n fc n I (l-n i )(l-n i )(l/2-2S fc .S,)J (1.4) 

This terms " counts" the number of singlet pairs of next- 
nearest neighbours (parallel dimer s) in all uncrossed pla- 
quettes. The resulting Hamiltonian Hw + Hk has a 
structure similar to that of the Roksar-Kivclson QDM, 
with both terms flipping dimers and terms counting the 
flippablc pairs of dimers. Hw can also be interpreted as 
a 4-site ring-exchange term on uncrossed plaquettes [l(| • 
It also presents a RK point (here at W = t 2 ), while for 
W/I^l — ► — oo an ordering in chains is favored; varying 
Wjti allows to make a continuous connection between 
both these limits and the case W = 0, and to understand 
better the physics around this point. 

C. Purpose of the study: phase diagram 

Since t 2 and W are the only energy scales in this model 
(at zero temperature) we aim at determining the phase 
diagram as a function of the ratio W ' jt 2 (or W, if we set 
t-2 = 1). First, we notice that for t 2 = W, the Hamil- 
tonian has the same property as in the quantum dimer 



model at the RK point: it can be written as a sum of pro- 
jectors (one per uncrossed plaquette). At this point (the 
RK point of the t 2 — W model) the wave function with 
an equal amplitude on all configurations (in each sector 
of connected configurations) is annihilated by each pro- 
jector, and thus is a ground state with zero energy. For 
W > t2, again for similar reasons as in the QDM, configu- 
rations of minimal energy arc those without any flippable 
pair of spins, and these configurations are (degenerate) 
ground states with zero energy. The center of interest of 
this study is the case where W < t 2 , i.e. the region be- 
tween the RK point and the W = — oo point, where the 
ground state is composed of decoupled Heisenberg chains 
(its energy is the sum of the energy of these chains with 
an AF coupling J = 2\W\ and a charge term L tot W/2 
where L to t is the cumulative length of the chains). In- 
stead of considering W and t 2 , one can define a reduced 
parameter 6 = BXctaax{W / 1 2 ) varying continuously be- 
tween the decoupled Heisenberg chains (6 = —tt/2) and 
the RK point (6 = tt/4). 

Between these points, the different phases one can ex- 
pect are (i) a columnar phase, ordered in chains, with 
rotational symmetry breaking, translational symmetry 
breaking in one direction (perpendicular to the chains) 
and thus a 4-fold degeneracy of the corresponding ground 
state; this phase is encountered for 9 = —ir/2, and could 
a priori extend over a finite range of 9 in the vicinity 
of that point. Note that the term columnar is used al- 
though this phase differs from the columnar phase of the 
QDM on the square lattice, due to the additional spin 
degrees of freedom; (ii) a RSPC or plaquette phase, with 
the full rotational symmetry of the lattice but a breaking 
of translational symmetry in both directions and a 4-fold 
degeneracy of the GS; (iii) mixed phases, with rotational 
and translational (in both directions) symmetry breaking 
(and a 8-fold degenerate GS), corresponding to a dimer- 
ization of the Heisenberg chains. A priori two types of 
mixed phases appear naturally depending whether neigh- 
boring chains dimerize in phase, or in antiphase. Note 
that these phases are a natural extension of the one re- 
cently discovered in the context of the RK QDM [l3| . We 
shall refer to them as Mixed Columnar- Plaquette Crystals 
(MCPC). 

Knowing that 9 = —tt/2 corresponds to a columnar 
phase, we refer in phase or in antiphase dimcrisation as 
MCPC- 1 and MCPC- 2, respectively, as shown schemati- 
cally in Fig. H 

These phases have distinct symmetries and the corre- 
sponding ground states are characterized by different sets 
of (4 or 8) quantum numbers, which we define using the 
following conventions: the x and y axes of the lattice are 
parallel to the links of the square lattice (on which the 
dimers live), and the unit length corresponds to one link 
of this square lattice. To define point group symmetries 
(those of the C^ v point group - or C 2v or C v for certain 
wave vectors) we set the center O of the lattice at the 
center of an uncrossed plaquette. The quantum numbers 
of the degenerate GS of the various phases are listed in 
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FIG. 2: (Color online) Possible scenarios for phase diagrams 
at n — 1/4 as a function of the ratio of parameters W/t?, 
depending whether an infinitesimal ti coupling dimerizes an- 
tiferromagnetic (AF) chains from a columnar order into a 
Mixed Columnar- Plaquette Crystal of type 1 (MCPC-l, up) 
or 2 (MCPC-2, down). The transitions indicated between ei- 
ther of these phases and a RSPC should be presumably of 
first order in the first case and second order in the second 
case, but one cannot exclude that the MCPC phase extends 
all the way to the RK point. 



D. Methods 



In the present work, we first discuss the regime near 
W/t2 — — oo perturbativcly (Section |TT|. Next, imple- 
ment a simple variational approach (discussed further 
in Section adapted to describing the various can- 
didate phases; in a second step, we shall use Lanczos 
Exact Diagonalization techniques to study the t<i — W 
model on clusters with periodic boundary conditions in 
both directions. The sizes of the clusters we consider in 
numerics are TV = 32,72 (7r/4-tiltcd checkerboard clus- 
ters corresponding to untilted square clusters of lengths 
L = 4, 6), TV = 48 (7r/4- tilted checkerboard cluster cor- 
responding to an untilted rectangular cluster of dimen- 
sions (L x ,Ly) = (4,6)) and TV = 64 (untilted checker- 
board cluster corresponding to a 7r/4-tilted square lattice 
of length L = Ay/2). Except for TV = 72, the cluster peri- 
odicity is compatible with all wave vectors q mentioned 
in the table above (for the TV = 72 cluster, the wave vec- 
tor q = (tt/2, 7r)-and those equivalent to it up to point 
group symmetries- are unaccessiblc). 



Phase — > 


MCPC-l 


MCPC-2 


RSPC 


Columnar 


(Al,9 = (0,0)) 


X 


X 


X 


X 


(Bl,q= (0,0)) 


X 


X 





X 


(Al,q= (7T,7T)) 


X 





X 





(A'l,?=(ir,0))(l) 


X 


X 


X 


X 


(Bl, g =(jr,jr)) 


X 











(A'l,g=(7r,0))(2) 


X 











(B',q = (tv/2,tv)) 





X 









TABLE I: Quantum numbers of the degenerate GS (in the 
thermodynamic limit) associated with the various phases ex- 
pected in the t-x — W model at n = 1/4 on the checkerboard 
lattice (for W < t^). The X sign indicates that a wave func- 
tion with the corresponding quantum number belongs to the 
degenerate GS manifold for the corresponding phase (0 other- 
wise). Irreducible representations labeled with ' are respective 
to subgroups C2v (A'l) and C v (A 1 , B') of the point group C^ v 
when the wave vector considered is non-invariant under d v . 



Table H 

Given the symmetries of the candidate phases, one can 
make a guess about the nature of phase transitions in 
the model, for both sccnarii described in Fig. O In the 
case of a MCPC-2 phase for |t 2 | <C — W, giving way to 
a RSPC close to the RK point, the transition between 
those should be of first order, since the symmetry groups 
of one phase is not included in that of the other. By con- 
trast, since the MCPC-l phase distinguishes itself from 
the RSPC by the breaking of one of its symmetries (in- 
variance by a 7r/2 rotation), the transition between those 
could be of second order. 



Note that we restrict ourselves to the sector S z = 
(which includes all total spin sectors) - for convenience. 
In addition, we consider only configurations for which 
the z-components of the total spin on each row (of ver- 
tical bonds) and each column (of horizontal bonds) of 
the square lattice (which are conserved quantities in the 
present model) are zero. This condition is satisfied by 
the ground states corresponding to any of the expected 
phases, and allows us to reduce the number of colorings of 
any dimer (=charge) configuration (hence the total size 
of the Hilbert space). By using the character of spin in- 
version, all point group symmetries and translations (in 
fact, due to the numerical technique for encoding config- 
urations, we use translations not interchanging the sub- 
lattices of the dimer lattice, hence TV/4 translations in- 
stead of TV/2) the number of representatives for N = 64 
and TV = 72 clusters are close to 4.5 • 10 4 and 1.9 • 10 5 
respectively [l6| . 

A powerful tool to determine the phase diagram is the 
analysis of the lowest energy levels in symmetry sectors 
associated to each of the quantum numbers mentioned 
(and with a character 1 for the spin inversion Sf — > —Sf). 
For the various phases one expects a degeneracy between 
the quantum numbers marked by X in the table HI in the 
thermodynamic limit (TV — > oo). On finite clusters, this 
degeneracy is lifted and the lowest state is always found 
for (Al, q = (0, 0)) - hence we look at the lowest excited 
states. Ideally, an unambiguous signal of spontaneous 
symmetry breaking is provided by the collapse of the cor- 
responding excitation energies AE = E t — E (A1, (0,0)) 
with increasing TV. However, the low-energy spectrum 
on large enough a cluster (e.g. TV = 64) gives enough 
information for a first analysis. 



II. PERTURB ATIVE APPROACH FROM 
W = -oo: COUPLING OF HEISENBERG CHAINS 

Adding the potential term W, to the Hamiltonian has 
several benefits, among them the existence of two partic- 
ularly simple special points, namely the RK point, and 
the point at W = — oo. For the bosonic model, the lat- 
ter yields the simple columnar configurations as ground 
states. However, W — — oo does not always present such 
a simple setting. For the triangular RK model, there are 
two families of ground states, each comprising a num- 
ber of members exponential in the linear system size, the 
degeneracy between which is not lifted until perturba- 
tions to leading non-trivial order in ti/W are taken into 
account 

The situation for our model is different still. At 
W = — oo, the ground state is obtained by maximis- 
ing the number of plaquettes with a pair of dimers in a 
singlet configuration. This leads to formation of a state 
breaking rotational and translational symmetries — just 
like the columnar state — but in which the spin correla- 
tions along the columns are critical. Indeed, W = — oo 
corresponds to decoupled Hciscnberg cains. 

The question appropriate to the setting of small \ti/W\ 
is thus: what is the most relevant perturbation induced 
by the kinetic term. This question has been addressed 
in — formally related — contexts by Essler, Tsvelik and 
Starykh + coworkers J3, EH ■ We closely follow the ap- 
proach of the latter [lj]. Their observation that in the 
Heisenberg chain not only the staggered spin but also the 
staggered energy correlations are critical — both decay 
as 1/r — is central: the chains are close to not only Neel 
but also to dimer ordering. 

For our model, the coupling of the staggered dimerisa- 
tion between neighbouring chains is symmetry allowed, 
and hence will generically appear as a perturbation is 
added. Indeed, it is easy to see how this happens in 
our model. For finite \t2/W\, flipping the dimers in two 
neighbouring plaquettes in adjacent rows yields an en- 
ergy gain of 0(W) for the plaquette marked by a cross 
in Fig. 02 whereas there is no such gain for the two plaque- 
ttes marked by circles, the coupling between the chains 
is thus generated at 0(|t2/W] 4 ), as each plaquette needs 
to be flipped out of the chain and back. 
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As analysed in Ref. [l4 |. this coupling is relevant and 
it will immediately lead to an in-phasc dimerisation of 
adjacent chains. In our above classification, this corre- 
sponds to a MCPC-1 phase, which is thus present in the 
limit W/t 2 — > — oo. 

Before we move on to variational (Section IIII[) and nu- 
merical ( Sections II VI andlV |) investigations of the effective 
t 2 — W model on the checkerboard lattice, let us make 
some further remarks on the regime where Wjti <C — 1. 
This will serve as a warm-up exercise for identifying 
ground states by quantum numbers. Let us consider 
a pair of Heisenberg chains of length L (with periodic 
boundary conditions, the chain axis being x), of spins- 
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FIG. 3: (Color online) Origin of dimerisation in the limit of 
large negative W (singlets flipped by Hk are represented by 
an up-down spin pair): the plaquette marked by a cross gains 
an energy W/2, while there is no such gain on plaquettes 
marked by circles. 



-^Si— <£^3— ^Lis>-- 



K/J 



i-i-hi-l-i- (a.f.) 

FIG. 4: Types of dimerization occurring in 2 Heisenberg 
chains coupled by a 4-spin coupling: in-phase dimerization 
(K < 0) or antiphase dimerization (K > 0) 



1/2 S^i and Si. 2- The Hamiltonian, 
H = J J2 S + - K "£(Si,i-S 

i+i,i)(Si ) 2-Sj + i ) : 



captures then coupling described above, although the 
true effective coupling is more complicated, in particu- 
lar, it includes similar four-spin terms along the chain 
favouring the same dimerisation pattern. 

At K = 0, the ground state is the product of ground 
states of the Heisenberg model on each chain, and has a 
symmetry (A, k x = 0)- A (B) labeling even (odd) states 
w.r.t. the chain interchange. The first excited state in 
the 5=0 sector correspond to a 2-triplet excitation on 
one chain (the other chain remaining in the GS). It is 
doubly degenerate, the quantum numbers of both states 
being (A,k x = 7r) and (B,k x = n). This state has (sim- 
ilarly to 1-triplet excitations) an excitation energy pro- 
portional to 1/L, thus collapsing to the ground state in 
the thermodynamic limit. 

When a weak interchain coupling K is added, the de- 
generacy of this first excited state is lifted, with a split- 
ting proportional to K/J, as shown on Fig.0 Since each 
of these states has the symmetry of a dimerized state, 
with dimerization either in phase (state (A,n)) or in an- 
tiphase (state (B,ir)), the sign of this splitting is there- 
fore associated with the type of dimerization susceptible 
to spontaneously appear in the system: for K/J > 0, 
the (B, 7r) state has a lower energy than the (A, tt) state 
and the system tends to dimerize in antiphase; this is the 
opposite case for K/ J < 0. 
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Let us examine now the finite size scaling of these ex- 
citations in presence of the 4-spin coupling. The excita- 
tion energy of the lowest excited singlet state ((A,ir) for 
K < and (B,ir) for K > 0) vanishes in the thermo- 
dynamic limit (the convergence is as 1/L for K/J C 1), 
while the excitation energy to the lowest triplet converges 
to a finite value (spin gap) (see right plot on Fig. [6J. Con- 
sequently, while on a small enough system the triplet ex- 
citation has lower energy than the lowest singlet, these 
levels cross as a function of system size (at fixed K/J). 
This is illustrated (in a slightly different way) by the 
left plot on Fig. [HI showing the positions of crossings 
(K/J) + /_(L) between the lowest triplet and either the 
lowest (A, 7r) or (B,ir) singlet, as a function of system 
size. One can see that these values converge to zero 
as L — > oo. (Rigorously, one cannot be fully conclu- 
sive with the present data set but a vanishing value as 
L — > oo is expected from the reasoning above. Moreover 
we checked that expressions of the type C/L + C / L 2 
fit the data of (K/J) + /_ better than any expression of 
the type C + C /L a .) Therefore the lowest singlet ex- 
citation, of symmetry either (A, ir) or (B,tt) depending 
of the sign of K, collapses onto the GS in the thermo- 
dynamic limit while a spin gap survives above. This is 
the well-known scenario of a spontaneous dimerization, 
which is precisely the type of scenario one expects in our 
two-dimensional effective model in the limit of weakly 
coupled chains (of colored dimers). Here, it is important 
to note that there are distinct reasons for the vanish- 
ing gaps. Even for W = —oo, there will be finite-size 
gaps of 0(|1F|), which vanishes algebraically due to the 
criticality of the chains, while the other gap, paramet- 
rically small in \t 2 /W\, collapses due to the presence of 
symmetry-breaking . 



III. VARIATIONAL APPROACH DESCRIBING 
THE CANDIDATE PHASES 



A. Principle and trial wave functions 
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FIG. 5: (Color online) Energies of the ground state and first 
excited states of two periodic Heisenberg chains of L = 10 
spins 1/2 (Heisenberg coupling J) coupled with a 4-spin in- 
terchain coupling K. (K/J) c indicates the crossing between 
the lowest triplet and the (A, k x = tt) singlet. Notice also 
the degeneracy of (B,k x = it) and (A,k x — 0) states at the 
Majumdar- Ghosh point K/J — 4/3. 
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Before analyzing the exact ground state and lowest ex- 
citations on finite systems, we estimate the energies of 
trial wave functions associated with the different candi- 
date (RSPC, MCPC-1 and MCPC-2) phases in the range 
of parameters of interest (— ir/2 < 8 < tt/4). A compari- 
son between their variational energies provides informa- 
tion on their relative stability. The trial wave functions 
we consider are built as the tensor product of all equiva- 
lent plaquettes of an identical wave function \ip p ) defined 
on a single uncrossed plaquette p, with a resonating sin- 
glet dclocalized on the 4 sites of the plaquette. 

On this plaquette, \if} p ) is expanded over the 4 S z = 
configurations with 2 particles on the plaquette (and 
respecting the dimer constraint): if the 4 sites of the 
plaquette are labeled from 1 to 4 clockwise around the 



FIG. 6: (a) (Color online) Position of the energy crossings 
between the lowest triplet and lowest (A, it) ((K/J)-) and 
(B,tt) ((K/J) + ) singlet states as a function of 1/L (L even 
up to 16); (b) Excitation energy AE of the lowest triplet for 
different values of K/J versus 1/L. Dashed lines are only a 
guide to the eyes. 



plaquette (starting from e.g. the site on the upper left 
side) these 4 configurations are: \u) = | fi O2 I3 O4) 
(0 denoting an empty site); \d) = \ J4 O2 T3 O4); 
\r) = |0i T2 3 U) and |/) = |0i | 2 3 T 4 ). M<t>) « a 
linear combination of the 2 sing let states (\u) + \d))/V2 
and (\r) + \l})/V2 (the plus sign results from the unitary 
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transformation mentioned in II B[) : 

i , \ / <\ ,,J«) + \d) . ,,s \r) + |0 
= cos(0) ' 71 7 + sin(^) 1 71 7 

. The parameter ^ can be restricted to values between 
and 7r/4, thus describing mixed phases obtained by 
dimerization of horizontally oriented Hcisenberg chains. 
The global wave function on a iV-site cluster \^>o(<p)) = 

( 3i'Ji\'ipp i {<fi)) depends on the angle <f> and the set of 
N/8 chosen plaquettes pi, i.e. the type of dimerization 
described: either in phase (MCPC-1) or in antiphase 
(MCPC-2). In the first case the RSPC corresponds to 
an angle <fi = it/ A. The wave function |^o) need not be 
symmetrized w.r.t. space group symmetries, in order to 
compute of the expectation values of Hk and Hw in the 
thermodynamic limit: indeed, for two symmetry-related 
and distinct vectors l^o) and I ^o), quantities such as 
(^' \H k i w \^q) give a relative contribution 0(1/N). 



B. Trial function for MCPC-1 and RSPC phases 

In order to describe the MCPC-1 phase and its sta- 
bility w.r.t. the RSPC, the trial wave function |>]/J)(</>) 
is used to compute expectation values of Hk and Hw 
The only terms of Hk contributing to (^\\Hk\^o) are 
plaquette flips on occupied plaquettes, exchanging u/d 
configurations and r/l ones. The average kinetic energy 
(per particle) is then: 

H K ,i{<t>) = ±/N{*l\H K \*l) 

= -2t 2 sin(» cos(0) (3.1) 

For the non-diagonal part of Hw (with Sf SJ terms), 
again only terms for which both sites are in an occupied 
plaquette contribute. But concerning the diagonal part 
of Hw, terms 1/2 — 2Sf 5? with i, j on a void plaquette 
between 2 occupied plaquettes also contribute, propor- 
tionally to either cos(</>) 4 or sin((/>) 4 depending on the 
position of this void plaquette. The expectation value of 
Hw (still per particle) is: 

H w , 1 {ct>) = A/N{^l\H w \^l) 
= W(l + sin(0) 4 /4 + cos(0) 4 /4) (3.2) 

The minimization of (Hk(4>)) + (Hw(4 > )) w.r.t. </> can 
lead to 2 distinct cases: 

For W /ti < —4 a non-trivial value 4>\ of </> minimizing 
the expectation value of H = Hk + Hw is found, corre- 
sponding to a MCPC-1 state breaking the 7r/2-rotational 
symmetry. <fii is solution of: W/t2 = — 4/ sin(</>i), and 
the corresponding average energy is: 

E 1 ( ( jy 1 ) = 2^ + 5W/4 

For W/t 2 > —4 the minimization gives cf>i ~ 7r/4, 
which means that the rotationally invariant RSPC is the 



most favorable state in this approach. The average en- 
ergy of the RSPC is estimated as: 

£i(7r/4) = -t 2 +9W/8 

The energies Ei(n/4) and (for 9 < arctan(-4)) Ex{<j>i) 
are shown as a function of the parameter 9 in figure [7] 
Notice that for 9 < arctan(— 4), Hi(tt/A) corresponds 
to a local maximum of the function Hi (</>) and is shown 
only for comparison to the (physically relevant) MCPC-1 
variational energy E\ (</>i ) . 




FIG. 7: (Color online) Energies (per particle and in units 
of \Jt\ + W 2 ) of variational ground states describing the 
MCPC-1 [9 < arctan(-4)], RSPC and MCPC-2 phases as a 
function of 9. Exact ground state energies for N = 32, 48, 64 
and 72 clusters are shown for comparison. 



These results give also an indication of the nature of 
the MCPC-1 /RSPC transition: the expansion of the 
variational energy functional Hk,x + P^w,i{4>) m powers 
of x = (f> — 7r/4, for W/t2 close to —4, is: 

Hx(4>) = i7 1 (7r/4) + 2(* 2 + IF/4).T 2 + (-19IF/24-2t 2 /3)x 4 

The sign of the coefficient of the x 2 term changes for 
W = —4*2 while that of the x 4 term remains positive 
around that point: in the frame of Landau's theory of 
phase transitions, this is characteristic of a second order 
transition with x as an order parameter, which varies 
continuously around the MCPC-1/RSPC transition and 
vanishes in the 7r/2-rotationally invariant phase. 



C. Trial function for the MCPC-2 phase 

In this case, the trial wave function \^>q) is still param- 
eterized by an angle 4> between and 7r/4 and differs only 
from the wave function describing the MCPC-1 phase by 
the position of the occupied plaquettes. The angle 7r/4 
corresponds here to a rotationally non-invariant pattern 
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of rotationally invariant occupied plaqucttcs (which dif- 
fers from the RSPC). Again, the structure of the trial 
state implies that the contribution to the expectation val- 
ues of Hk and the non-diagonal part of Hyy comes only 
from terms acting independently on occupied plaquettes. 
But here the diagonal part of Hw acts differently on 
the wave function than in the previous case: it gives a 
non-zero expectation value only on the void plaquettes 
situated between 2 occupied plaquettes to the right and 
left (see Fig. [2]). The corresponding term in (H w ) is 
thus proportional to the probability cos(0) 4 for particles 
of both plaquettes to be in \r) or \l) states. Eventually 
the expectation value of H as a function of </> reads: 



H 2 (<j>)=4/N{Vl\H w + H K \*l) 
W(l + cos(0) 4 /4) - 2t 2 cos(</>) sin(<£) 



(3.3) 



and is minimized, either for (f>2 = 7r/4 when W > 0, or 
when W < for 62 solution of: 



W/t 2 = -4 



1 + tan( 



tan(202 



The angle </>2 (for W < 0), and consequently the corre- 
sponding expectation value £2(^2) of H , has no simple 
expression as a function of W jti or 9; a numerical reso- 
lution proves that for W < it is greater than Hi(cf)i) 
found with the wave function (see Fig. [7|) - this is 
not a surprise since the (W/4) sin((/>) 4 term present in 
Eq. 13.21 is absent in Eq. 13.31 In other words, this ap- 
proach indicates that the MCPC-1 or RSPC phase is 
stabilized w.r.t. the MCPC-2 phase by interactions on 
some plaquettes, as soon as they are attractive (W < 0). 

For W > the variational ground state |^q) corre- 
sponding to an angle $2 = 7r/4 has an energy 17VK/16 — <2 
lower than that of the RSPC found before. Hence this 
can be considered as the variational ground state in this 
approach, predicting a domain of stability — 4t 2 < W < 
for the MCPC-1 phase (see Fig. [8]) - but one has to take 
into account the limitations of this approach, discussed 

innm 



D. Comparison to the bosonic case 

In this paragraph we apply the previous variational 
method to the bosonic case, i.e. to the QDM of Roksar 
and Kivclson on the square lattice, to have a comparison 
between the variational and exact phase diagrams. In 
this variational wave function \^)(<f)) describing 

the plaquette phase and the mixed phase (both described 
in [13], i.e. the bosonic analog of the MCPC-1 phase ; 
similarly a mixed 2 phase can be defined as the bosonic 
analog of the MCPC-2 phase) is defined, as previously, as 
a product of local wave functions on plaqucttcs occupied 
in a plaquette pattern: 



(\v) and \h) correspond to either 2 vertical or 2 horizontal 
dimers on the given plaquette). The energy per particle 
of the state l^)^) is here (—t and V being the amplitudes 
of the kinetic and potential terms in the QDM, defined 
as in [IH): 

Hi(<j>) = -2£sin(^)cos(0) + V(l + cos 4 (0) + sin 4 (0)) 

For V > — t this is minimized for <j) — 7r/4 which corre- 
sponds to a plaquette state, while for V < — t the angle 
(f>i is such that V/t = — 1 / ' sin(2<f>i) . The corresponding 
energy per particle is: 



t- 



V 



AV 

£i(7r/4) = 3V/4-4/2 



(V/t<-l) 

(v/t > -1) 



|* p ) = cos(<j>)\v) + sm(<f>)\h) 



As in the fcrmionic case, one can define a variational wave 
function paramctered by an angle <f>' and corresponding 
to the mixed 2 phase, and show that, after minimization 
w.r.t cf) 1 it has a lower energy than £"i(7r/4); but again, 
in that case the variational method is not adapted to the 
situation close to the RK point. 

Consequently, the variational approach for bosons on 
the checkerboard at n — 1/4 predicts the existence of a 
mixed phase for V < —t and a transition (second order 
again) to a plaquette phase at V = —t. In this approach 
the plaquette phase extends from that point up to V = 0, 
while the ground state energy in the domain < V < t 
is better approximated with a mixed U-type variational 
wave function- this is essentially due to the inadequacy 
of this method close to the RK point. The variational 
phase diagrams for both the fermionic and the bosonic 
model (both for an average occupation number n = 1/4) 
are shown on Fig. [8] 

The larger extent of the RSPC phase in the fermionic 
model than of the plaquette phase in the bosonic can 
be attributed to the fact that, in the first case, a sin- 
glet resonating on a plaquette in a RSPC phase is less 
coupled to neighbouring occupied plaquettes: the pro- 
cess coupling neighbouring plaquettes is subject to the 
constraint that interacting particles have opposite spin, 
this constraint being absent in the bosonic case. Con- 
sequently, for fermions the RSPC phase is more stable 
and the transition to a MCPC phase with longer range 
correlations occurs for a larger (negative) value of W. 



E. Reliability of the variational approach - 
comparison to exact ground state energies 

Let us now comment briefly on the reliability of this 
variational approach to give a qualitative, or even quanti- 
tative, estimation of the phase diagram. For this, we have 
compared the variational ground state energies Hi(<pi) 
and H2(<fi 2 ) to the exact ground state energies obtained 
on periodic clusters of size N = 32,48,64 and 72 (see 
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W=-4£/ 



(0=-7t/2) 




RK(W=t 2 ) 




V=0 

FIG. 8: (Color online) Variational phase diagrams for the 
model of fermions (up) or bosons (down, with notations of 
ref. [l3|]). The point of transition between a 7r/2-rotational 
invariant phase and a phase breaking that symmetry is pre- 
dicted at W = — 12 (for fermions) and V — — t (for bosons). 
For < W < f 2 (resp. < V < t) the MCPC-2 or mixed-2 
phase is the most competitive of our crystalline trial wave- 
functions, but it in turn loses out to the simple RK wave 
function, see IIII El (green vertical lines). 



Fig. [7J. Since the orientation and geometrical shape 
of these clusters differ from each other, one cannot do 
an accurate finite size scaling of the exact ground state 
energy that would allow for a precise comparison be- 
tween exact and variational energies in the thermody- 
namic limit. However, it appears clearly that finite-size 
effects on the exact energies (per particle) decrease with 
N, which suggests that for a wide range of 9 (between 
— 37r/8 and 0) the exact and variational ground state en- 
ergies differ from about 5% or less, and their variations 
with 9 are very similar [l7| . The discrepancy between 
exact and variational results is larger (?) when 9 gets 
close to 7r/4: at the RK point, the ground state energy 
on any cluster considered is zero while the variational 
energy per spin is £2/16 and corresponds to a MCPC- 
2 phase unexpected here: the failure of the variational 
approach close to the RK point originates in the inad- 
equacy of the variational wave function when the exact 
ground state is much more disordered (which holds also 
in the bosonic case); there, the RK wave function is a 
better trial wave function, with an energy per particle 
E RK =W-t 2 [(V - f)/2 in the bosonic case] [H|. This 
RK trial wave function has an energy lower than |^q) 
for either W/t 2 > or V/t > (see Fig. [8j), defining a 
domain of stability of a RK phase, which has to be inter- 
preted more appropriately as a domain where the varia- 
tional approach fails; (ii) in the 1 2 — W model when 9 is 



close to — it/2, the failure of this variational approach is 
expected since the exact ground state consists of Heisen- 
berg chains (weakly dimerized for finite t 2 ), with spin- 
spin correlations along chains slowly decaying and the 
trial wave functions describing isolated resonating pla- 
quettes are no more valid there. For that reason, the 
position of the MCPC-1/RSPC transition, predicted in 
this approach at W/t 2 = —4 (hence 9 = —1.326(1) rela- 
tively close to — tt/2) can differ appreciably from the real 
position of this transition. To determine the latter, one 
has to treat the model exactly, taking all allowed config- 
urations into account (and not only those characteristic 
of the plaquette ordering) by methods such as exact di- 
agonalisation. 



IV. ANALYSIS OF THE LOW-ENERGY 
SPECTRUM AND IDENTIFICATION OF 
PHASES 



In the limit t 2 <C \W\ of the t 2 — W model, where the 
variational approach of Sect ion [ITT] is least reliable, the ki- 
netic term Hk can be described as a perturbation, while 
the unperturbed Hamiltonian H\y has for ground state 
a product of Heisenberg chains. The effect of Hk can be 
described by an effective coupling K e g W\t 2 /W\ 4 , as- 
sociated with processes of order 4 out of the Heisenberg 
ground state, where 2 singlets belonging to 2 neighbour- 
ing chains and opposite to each other are flipped and then 
flipped back; the interaction resulting from this process 
is attractive due to the Hw term on the interchain pla- 
quette^), making this process the most important one in 
perturbation in t 2 /W (among those having an influence 
on the type of dimerization). Although we did not de- 
termine analytically the sign of this effective coupling for 
< t 2 <C \W\, by analyzing the energy splitting of the 
the analogs on the checkerboard of the (A, it) and (B, tt) 
excitations (of the 2-chains J — K model) we determine 
the sign of K c g, i.e. the type of dimerization occurring 
in the system. 



Weak-coupling regime: low-energy spectrum 
and quantum numbers 



Let us consider first the t 2 = limit: here the low- 
energy spectrum has a simple structure, i.e. lowest states 
are composed of lowest states of Heisenberg chain of the 
corresponding length (3 chains of L = 6 for the N = 72 
cluster, 2 chains of L = A.k for the N — 32. k cluster (k = 
1,2). The ground state is found in sectors (A1,(0,0)), 
(51,(0,0)), (A'1,(tt,0)) and (A'l, (0, tt)) (the 4-fold de- 
generacy accounts for the 4 ways of regularly accommo- 
dating 2 Heisenberg chains on the checkerboard cluster). 
At t 2 = 0, and by extension in the weak-coupling regime 
\t 2 /W\ <§; 1, the first excited states are, either on the 
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FIG. 9: (Color online) Energies of lowest excited states (in 
the S = sector) of symmetries Al/Bl,q = (0, 0)/(n, 7r), 
(A'l, (tt, 0)), (B 1 , (tt/2, tt)), at ti < | W|, for a cluster of N = 
64 sites. Inset: splitting between a 2-triplet (B' , (-k /2,n)) 
state (plus symbol) and other 2-triplet states (X symbols). 



cluster N = 64 (see Fig. %B or N = 72: (%) a state in 
the 5=1 sector (degenerate between various quantum 
numbers) corresponding to a 1-triplct excitation on one 
chain and the Hciscnbcrg ground state of the othcr(s) 
chain(s); (ii) a state in the sector 5 = 0, with all quan- 
tum numbers listed in the table HI corresponding to a 
charge excitation breaking of a Heisenberg chain into an 
isolated singlet and an open (L — 2) chain); (Hi) a state 
in the sector 5 — corresponding to 2 triplet excita- 
tion on the same Heisenberg chain, the other chains be- 
ing in their ground states. On the N = 64 cluster this 
state is found with the quantum numbers (Al, (jr, 7r)), 
(Bl,(7r,7r)), (A'l,(7r,0)), (B',(tt/2,tt)) (and those re- 
lated by symmetry); (iv) slightly above the latter, a state 
in the sector 5 = corresponding to two 1-triplet exci- 
tations on distinct chains. 

Although the lowest excitation in the 5 = sector 
is, on clusters considered, the charge excitation, as far 
as these states can be labeled as charge and 2-triplet- 
excitations we rather focus on 2-triplet excitations (states 
(Hi) and (iv)), since their excitation energy is a finite-size 
effect, being proportional to \ W\/L which vanishes in the 
thermodynamic limit, while state (ii) has an excitation 
energy of order \W\ even for L — > oo. Among 2-triplct 
excitations, state (Hi) (2 triplets on the same chain) is of 
greater interest: not only does it have a slightly lower en- 
ergy than state (iv), but a reasoning based on the model 
situation of 2 Heisenberg chains coupled by a 4-spin cou- 
pling term (see Sec. [TTJ) shows that this state gives infor- 
mation about the type of dimcrization favoured by the 
interchain coupling (or £2 here). 

To analyse the influence of a weak kinetic coupling 
£2 <C |W| on the system and determine, in the light of 
previously discussed features of the J — K model of cou- 



pled Heisenberg chains, we focus on states with 2 triplet 
excitations on the same chain. The states corresponding 
to this excitation with different quantum numbers split 
when £2 increases; at lowest order in t2/\W\ the split- 
ting occurs between the state in (B 1 , (n/2,ir)) (2D ana- 
log of the (B, tt) 2-triplet state in 2 coupled Heisenberg 
chains, hereafter B state) and states with other quan- 
tum numbers (analogs of the (A, n) 2-triplet excitation 
for 2 coupled Heisenberg chains; hereafter A states). In 
this splitting the B state (associated with a dimeriza- 
tion in antiphase) has higher energy than the A states 
(associated with a dimcrization in phase); moreover, we 
have checked that the energy difference between the B 
state and the A states scales as Wfo/W) 4 , as expected 
from the previously discussed comparison with the J — K 
model of 2 Heisenberg chains (where the corresponding 
splitting occurs linearly in Kj J) . Hence we conclude 
that on the checkerboard, for small but finite £2/|W| the 
Heisenberg chains dimerize in phase and form a MCPC-1 
phase, which is consistent with results of Section IIIIl 

While in the limit of small t2/W the physics of the 
model is quasi-lD, making finite-size effects very impor- 
tant on the checkerboard clusters analyzed (especially 
at £2 = 0, Heisenberg chains being critical), for larger 
couplings these finite-size effects become less relevant, as 
soon as the clusters can accommodate the various or- 
dered phases. Hence in the latter case and with the clus- 
ter sizes available, it is reasonable to analyze the t 2 — W 
model non-perturbatively, and without making reference 
to an effective model (such as that of coupled Heisenberg 
chains). 



B. Non-perturbative analysis for intermediate 
W/t 2 : MCPC-1 or RSPC phase ? 

Away from the £2 <SC \W\ limit, one cannot simply iden- 
tify each of the first excitations as 2-triplet or charge ex- 
citations, but by using symmetries of the model one can 
characterize them by their quantum numbers. As seen in 
Sec. II CI for each of the candidate phases, in the thermo- 
dynamic limit, we know the quantum numbers associated 
with wave functions of the (4- or 8- degenerate) ground 
state. Hence the relative order (in energy) of first exci- 
tations with those quantum numbers, for large enough 
systems, should be characteristic of the symmetry of the 
ground state, and thus of the phase in question. From 
Scc. lIV Al the MCPC-1 phase is expected to extend over 
a finite range of W '/ti, going either to the RK point or to 
a non-trivial point of transition towards a RSPC phase 
(in analogy with the QDM on the square lattice). In 
Fig. [10] we plot the lowest excitation energies (in units of 
y/t\ + W 2 ) corresponding to quantum numbers listed in 
Table |TJ as a function of the parameter 9 = arctan(VF/£2) 
(N.B. This is equivalent to considering the Hamiltonian 
H = sm(6)[H w (W = 1)] + cos(9)[H K (£ 2 = 1)] ). 
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FIG. 10: (Color online) Excitation energies AE = E / 1 — 
Egs (Ei if labeled by (*)) in units of + W 2 for different 
quantum numbers, in the singlet [l!| magnetic sector (open 
symbols) as a function of 9 = &rcta,n(W /t2), on the N — 
64(a) and N = 72(b) checkerboard clusters. The lowest triplet 
excitation is also shown (star symbols). 



The most striking feature of these graphs con- 
cerns the states with quantum numbers (Al, (tt, w)) 
and (51,(0,0)). The excitation energy (in units of 
\Jt\ + W 2 ) of the former, being that of the charge exci- 
tation in the 9 — ► —ir/2 limit, collapses when 9 increases. 
In comparison, the (SI, (0, 0)) state, degenerate with the 
ground state at t% — (uncoupled Heisenberg chains) 
becomes separated energetically from the ground state 
when t2 becomes non negligible (the vanishing of both 
excitation energies when W/t2 gets close to 1 is due to 
the degeneracy at the RK point). On the N = 64 cluster, 
for a wide range of 9, the Al, (n, n) state has an energy 
significantly lower than the (SI, (0, 0)) state and close to 
that of lowest (Al, (0, 0)) and (A'l, (tt, 0)) states (This is 
less obvious with data from the N = 72 cluster, where the 
effective length of chains L = 6 results in stronger finite- 



size effects). The eigenvalue crossing between (Al, (tt, tt)) 
and (51,(0,0)), similar to that observed for n = 1/2, 
indicates a breaking/restoration of rotational symmetry, 
signaling a transition between the MCPC-1 phase (ex- 
pected in the thermodynamic limit for t?. <C I14 7 !) and 
the RSPC for which the symmetrized wave functions 
of the ground state have quantum numbers (Al, (0,0)), 
(A'l, (tt,0)), (A'l,(0,7r)) and (SI, (tt, tt)). As the data 
of the N = 72 cluster are less clear, the present analysis 
should be extended to larger clusters, which is unrealistic 
with the current most advanced computational resources 
[2fj| unless we find other tools to be more conclusive about 
the existence and position of a MCPC-1/RSPC phase 
transition. 



V. PLAQUETTE CORRELATIONS ON THE 
CHECKERBOARD 

Since the finite-size effects encountered for clusters of 
sizes TV < 72, particularly in the weak-coupling regime 
and at the supposed MCPC-1 /RSPC transition, make it 
difficult to identify clearly this transition by analysing the 
low-energy spectrum only, one needs complementary in- 
formation about the nature of the ground state. The ED 
numerical technique employed in this study allows also 
to compute expectation values of observable and their 
associated correlations. Considering the structure of the 
different ordered phases expected in this model (both the 
MCPC-1 phase and the RSPC are composed of singlets 
localized on uncrossed plaquettes of the checkerboard), 
it is preferable to consider an observable defined on an 
uncrossed plaquettc (and then symmetrized) rather than 
on a single site. Hence in this section we discuss results 
about two types of plaquctte correlation functions in the 
ground state of the t^ — W model, computed on clusters 
N = 32,64,72. 

The first type of plaquctte operators for which we com- 
pute correlations have SI point group symmetry, and are 
related to the flippability of a given plaquette located at 
r (i.e. of an uncrossed plaquette of the checkerboard): 

P_(r) = (n lth n 2t h - ni tV n2, v )ds liZ+ s 2>z (5.1) 

The correlations of P_ operators are computed from the 
ground state wave function |^gs) (first determined as 
the ground state in the Al,q = (0,0) symmetry sector 
and then expanded on all configurations). The average 
value of P- on the ground state being zero by symmetry 
(for any plaquette), the correlation function is defined as: 

C-(R) = n R .(^ GS \P-(0)P-(r)\^ G s) 

where P-(0) is computed on a reference plaquette, P-(r) 
on a plaquette at distance R from the reference one and 
tir is the number of plaquettes at distance R from the 
reference. A non-zero value of C-(R) in the thermody- 
namic limit and at large distances corresponds to a phase 
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breaking the 7r/2-rotational symmetry (hence in this con- 
text to the MCPC-1 phase). 

In Fig. [TTJ are shown, for each of the cluster sizes con- 
sidered, correlations C a (R = d max ) at the maximal dis- 
tance between plaquettes equivalently occupied in a pla- 
quette phase (d max = 2\[2 for TV = 32, 72 and d max = 4 
for TV = 64). We also plot Fourier transforms S- of C_ 
correlations, for wave vectors q = (0,0) and q = (tt, n). 
These Fourier transforms are computed with a truncation 
at short distances, i.e. only correlations for distances 
r > 2 are taken into account (distances R = 0, 1 arc 
discarded since the plaquette operators on neighbouring 
plaquettes share at least one site; at distance R = \J~2 
the simultaneous double occupancy of both plaquettes is 
forbidden by the dimer constraint). 

The correlations of the Bl plaquette operator (P_) 
vary significantly with 9 between the Heisenberg- and 
the RK- limits, and give important information about 
the evolution of the ground state. Correlations of P_ de- 
crease strongly for 9 < — 37r/8 (which means W/t2 < 
W c ~ —2.4), both at the largest distance d max be- 
tween equivalent plaquettes, and in Fourier space at 
wave vector q = (0, 0) -in the weak-coupling regime (for 

9 < 9 C 1.2(1)) S-(q = (0,0)) is well approximated 

by a Gaussian function of 9 + tt/2; for 9 > 9 C , the ratio 
SL(0)/<SL(— 7r/2) becomes smaller when d max increases, 
so one can expect Bl plaquette correlations to vanish 
in the thermodynamic limit. This indicates that the ro- 
tational symmetry of the lattice, broken in the weak- 
coupling regime, is restored for 9 > 9 C , which is a signa- 
ture of a transition to a RSPC phase. 

We also computed (off-diagonal) plaquette-exchangc 
correlations, i.e. correlations of the kinetic operator 

Po(r) = (b\ 1 b] L + fe l t i 6] T )(& fcT 6 u + b kl b n ) + c.c. 

(where sites i,j,k,l are those around the void plaque- 
tte at position r). The connected correlations C p i(r) = 
(Pn(r)fb(O)) - (Pn(0)) 2 are vanishing in the 9 -> -tt/2 
limit, where charge moves away from Heisenberg chains 
are energetically forbidden; at the contrary, in a RSPC, 
they are expected to be important between resonating 
plaquettes (at relative position (2.p,2.q) — p,q 6 Z from 
each other) and significantly smaller otherwise. Con- 
sequently we focus specifically on plaquette-exchangc 
structure factors S p z(q) (Fourier transforms of the cor- 
responding correlations) at q = (%, 0), q = (0,vr) and 
q = (tt, 7r). They are shown for the TV = 64 cluster, 
along with the correlation in real space at maximal dis- 
tance C p i(d max ), on Fig. [TT] (graph (b)). 

These correlations increase significantly with 9 in 
the weak-coupling regime, indicating the appearance of 
resonating plaquettes characterizing the MCPC-1 and 
RSPC phases. We have checked (not shown) that the 
(7r, 0) structure factor is very well approximated by its 
contribution from correlations between resonating pla- 
quettes only, which indicates that the picture of a plaque- 
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FIG. 11: (Color online) Correlations in real space (C_) at 
distance d max (dashed and dotted curves), and in reciprocal 
space (S-) at q = (square symbols) of Bl plaquette opera- 
tors, on clusters of size N = 32 (a), TV = 64 (b) and TV = 72 
(c), as a function of 6 = a,xct&n{W j 't-z) . On graph (b) are also 
shown, for the TV = 64 cluster, plaquette-exchange correla- 
tions C p i at distance d max and the corresponding structure 
factor S p i at q = (tt, 0) and q = (tt, n). 
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tte pattern describes well the ground state wave function 
at 9 > 9 C , where these correlations are important. The 
smaller value of the structure factor at q = (it, ir) com- 
pared to that at q = (tt, 0) originates from non-negligible 
correlations at short distances such as d = \/5. The 
comparison between the correlation C p i(d max ) rescaled 
by a factor N/ 8—1 (accounting for the number of pla- 
quettes included in the Fourier sum defining S p i and 
that should be occupied and resonating in a MCPC-1 or 
RSPC phase, and the structure factor S p i(ir,0) itself, is 
eloquent: it indicates that for 9 < 9 C , almost all contri- 
butions to the Fourier transform come from resonating 
plaquettes, and the correlations between these plaque- 
ttes are long ranged. For larger values of 9 (> 9 C ) the 
(negative) contributions from other plaquettes become 
non-negligible; but structure factors at both wave vec- 
tors remain significant, and indicate the robustness of 
the ordering on these plaquettes. Although plaquettc- 
exchange correlations alone are not sufficient to distin- 
guish the RSPC from the MCPC-1 phase, their change 
of behaviour around 9 = 9 Cl associated with features 
of the low-energy spectrum and Bl (diagonal) plaquette 
correlations, indicates the presence of a MCPC-1/RSPC 
phase transition. 



VI. CONCLUSION AND PERSPECTIVES 

In summary, we have considered S — 1/2 fermions 
on the checkerboard lattice, with an extended Hubbard 
model in a limit of infinite on-site and strong nearest- 
neighbour Coulomb repulsions, at the specific fractional 
1/8 filling. In this limit, constraints characteristic of 
the square lattice dimer model naturally emerge, with 
the links on the square lattice corresponding to the sites 
of the checkerboard. Moreover, the lowest-order kinetic 
process allowed in perturbation (from the infinite repul- 
sion limit) flips 2 particles around an uncrossed plaque- 
tte, recalling the kinetic term of the Quantum Dimer 
Model, and the analogy is reinforced by considering an 
extra term similar to the potential term of the QDM. 
However, here the spin degrees of freedom of particles 
play an essential role, since kinetic processes act only on 
singlet states on a plaquette. The model we have consid- 
ered makes a continuous connection between the physics 
of critical Heisenbcrg chains (occurring in place of the 
columnar limit of the QDM) and a RK-like critical point 
of the present model. 

Starting from a situation where the Heisenbcrg chains 
of the first case are weakly coupled, we have identified the 
leading order - in perturbation in t-^jW - term coupling 
neighbouring chains, and characterized it as a relevant 
perturbation for the Heisenberg chains, driving their in- 
phase dimerization in the thermodynamic limit as soon 
as the coupling is finite. The corresponding phase (Mixed 
Columnar- Plaquette Crystal -1 ) distinguishes itself from 
other candidate phases (RSPC and columnar) of the 



e=7t/4:RK 




W=0 



9=-7i/2 
ID chains 

FIG. 12: Phase diagram of the t 2 - W model of S = 1/2 
fermions at n — 1/4 (in the range of parameter —tt/2 < 
9 < tt/4 where 6 = - arctan(W/t 2 )). The MCPC-1/RSPC 
transition occurs at 9 C ~ — 3tv/8 ± 0.2 according to exact di- 
agonalization results (a value 9 C = arctan(— 4) is estimated 
variationally) . 



model by a lower symmetry, and its extent in parame- 
ter space is determined by a detailed analysis of the low- 
energy spectrum obtained by exact diagonalization, tak- 
ing lattice- and time reversal- symmetries into account. 
This analysis indicates that, when going towards the RK- 
like point, the MCPC-1 phase persists up to a transition 
into a 7r/2-rotational invariant Resonating Singlet Pair 
Crystal, the analogue of the plaquette phase of the square 
lattice QDM; this transition is confirmed by the compu- 
tation of various types of correlations between uncrossed 
plaquettes of the checkerboard, characterizing the pla- 
quette ordering, and the symmetries of the ground state. 
The qualitative features of this phase diagram are also 
found by a variational approach, also indicating that the 
transition between the MCPC-1 phase and RSPC should 
be of second order. In particular, the system described by 
the extcndcd-Hubbard model with strong repulsions con- 
sidered first appears to be in a Resonating Singlet Pair 
Crystal similarly to the corresponding model at quarter 
filling [|. 

An open question is to know what happens at small but 
finite doping from the 1/8 filled case: a possibility is that 
the system remains in a crystalline phase confining the 
doping particles (either holes or electrons depending on 
the type of doping); alternatively, the RSPC could give 
way to a phase with either deconfined doping particles 
or the formation of bound Cooper pairs (as it happens 
in systems at small doping from 1/4 filling and 1/2 fill- 
ing 21() that would be an indication of a superconducting 
or supcrsolid phase. 
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